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1. Introduction 

This paper is devoted to propagation of singularities in linear hyperbolic par- 
tial differential equations with non-smooth coefficients in the framework of the 
Colombeau theory of generalized functions [2l |3] . Such equations involve a nonlin- 
ear interaction of coefficient singularities with those in the solution which emanate 
from the initial data. In the past years, this problem has been investigated inten- 
sively in a series of papers, in which the coefficients and the solution have been 
viewed as elements of the Colombeau algebra of generalized function. 

So far, microlocal elliptic regularity is well understood for general equations 
and systems in this framework [SIIHIIIH]- Further, generalized Fourier integral op- 
erators have been introduced in order to describe the wave front set of generalized 
solutions [3 ini [TUJ [TT] . These methods have allowed to describe the propagation 
of singularities in scalar first order hyperbolic equations with generalized function 
coefficients [5l [121 E] ■ As a new phenomenon, the Colombeau wave front set has 
been shown to possess a more refined structure than the corresponding distribu- 
tional wave front set [22]. For a survey of these methods and relations to various 
classical approaches we refer to jl5| . 

In the special situation of hyperbolic equations and systems with piecewise 
constant coefficients (with possibly jumps across smooth hypersurfaces), the prob- 
lem may be interpreted as a transmission problem and can be solved classically by 
a piecewise distributional solution. In case of systems or higher order equations, an 
incoming singularity will split at the interface into a refracted and a reflected wave. 
While in many cases it has been shown that the Colombeau solution is associated 
with the distributional solution of the transmission problem [131 HO] , the task of 
identifying the refracted and reflected singularity by means of the Colombeau wave 
front set - directly in the generalized solution without considering the associated 
distribution - has remained open. 

The purpose of this paper is to demonstrate - in a simple model problem 
- that the splitting of singularities at an interface can indeed be observed in the 
Colombeau solution. At the same time, the effect depends on the interplay of the 
scale of regularization of coefficients and data. We consider the Cauchy problem 
for the one-dimensional wave equation with a discontinuous coefficient depending 
on time 

d^u - cUfdlu = 0, t > 0, .T e R, 

, (1-1) 

u\t=o = uo, dtu\t=o = ui, X- e M. 

We will seek solutions in the Colombeau algebra ^^([0, oo) x K) of generalized 
functions, which will be defined in Section [2] below. The Colombeau algebras are 
constructed as factor algebras of nets of smooth functions, depending on a regu- 
larization parameter. They are differential algebras; the space of distributions can 
be imbedded as a subspace using regularization by convolution with a Friedrichs 
moUifier. The initial data as well as the coefficient will be taken from the algebra 
^^(R), which contains the space ^'(R) of distributions. 
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As a model situation, we consider the case when c(t) = cq + (ci — CQ)H{t — 1), 
uq ~ and ui ~ S, where co,ci > 0, cq 7^ ci, H is the Heaviside function and 
5 is the deha function. The singularity of the initial data propagates from the 
origin in the two characteristic directions until it hits the coefficient discontinuity 
at time t — 1. As the coefficient has a discontinuity in time, the singular support 
is expected to split at time t = 1 into a transmitted and a refracted part. In the 
classical setting, the piecewise distributional solution exhibits this splitting. We 
will demonstrate the occurrence of the splitting effect in the Colombcau setting as 
well, using the generalized singular support. However, in the Colombcau setting the 
effect depends on the scale in terms of the regularization parameter. If the initial 
data and the coefficient are regularized on the same scale, both the transmitted and 
the refracted ray will be shown to belong to the Colombcau singular support. If the 
coefficient is regularized by means of a slow scale mollifier, only the transmitted 
ray belongs to the Colombcau singular support. 

In order to be able to observe this effect, we have to prove a new existence 
result for the Cauchy problem which does not restrict the choice of scale 

for the regularization of the coefficient. The equation can be rewritten as a one- 
dimensional first-order hyperbolic system with discontinuous coefficients in the 
Colombcau algebra of generalized functions. This transformation will be used in 
the existence and regularity results. We also demonstrate that the Colombcau 
solution has the piecewise solution of the transmission problem as its associated 
distribution. That is, the nets defining the Colombcau solution converge to the 
piecewise distributional solution. This latter argument relies on energy estimates. 

The paper is organized as follows: we recall the definition and basic properties 
of the Colombcau algebra J# in Section[2] In Section[31 we show that problem 
is uniquely solvable in the Colombcau algebra $#([0,oo) x R) without restriction 
on the scale of regularization of the coefficient (Theorem 13. ip . The problem of 
propagation of singularities is addressed in Section HI (Theorem 14. 1[) . The question 
whether the regularity of the coefficient affects that of the solution is discussed in 
Section [5] (Theorem 15. ip . In Section [HI we show that the Colombcau generalized 
solutions to problem (jl.ip with arbitrary distributions as initial data admit the 
piecewise distributional solutions of the corresponding transmission problem as 
associated distributions (Theorem 16. 



2. The Colombeau theory of generaUzed functions 

We will employ the special Colombeau algebra denoted by Q'^ in |14j . which was 
called the simplified Colombeau algebra in [T]. However, here we will simply use 
the letter Q instead. Let us briefly recall the definition and basic properties of the 
algebra Q of generalized functions. For more details, see [14]. 

Let J7 be a non-empty open subset of M''. Let C°°(r2)^°'^l be the differen- 
tial algebra of all maps from the interval (0, 1] into C°°{ft). Thus each element of 
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C°°(r2)'°^^l is a family {u^)ee{o.i] of real valued smooth functions on 57. The subal- 
gcbra £Mi^) is defined by all elements {u^)ee{Q.i] 

of C°°(r2)(''^il with the property 
that, for all K il and a G Nq, there exists p > such that 

snp\d^u'{x)\^0{e-P) as e ; 0. (2.1) 

The ideal A/'(f7) is defined by all elements (M'^)gg(o,i] of C°°(r2)^°'^i with the property 
that, for all K (e fi, a e Ng and q>0, 

sup\d^u%x)\^0(e'^) ase;0. 
The algebra G{^1) of generalized functions is defined by the quotient space 

The Colombeau algebra on a closed half space [0, od) x M is constructed in a similar 
way. 

We use capital letters for elements of Q to distinguish generalized functions 
from distributions and denote by (u'^)eg(o,i] a representative oi U € Q. Then for 
any U, V E Q and a € Nq, we can define the partial derivative d"U to be the 
class of (9"M'^)eg(o,i] and the product UV to be the class of (u'^t;'^)eg(o,i]. Also, 
for any U = class of (^^(i, s))ee(o,i] £ G{[0,oo) x K), we can define its restriction 
C/|t=o G G{^) to the fine {t = 0} to be the class of (m^(0, a;))eg(o,i] . 

Remark 2.1. The algebra G{^) contains the space £'{^l) of compactly supported 
distributions. In fact, the map 

/ 1-^ class of if *Pe k2)ee(0,l] 

defines an imbedding of £' (fl) into Q{fl), where 

and p is a fixed element of 5(R'*) such that / p{x) dx = 1 and / x°'p{x) dx = 
for any a G Nq, |a| > 1. This can be extended in a unique way to an imbedding 
of the space 'D'{il) of distributions. Moreover, this imbedding turns C°°(f2) into a 
subalgebra of Q{il). 

Definition 2.2. A generalized function U £ Gi^) is said to be associated with a 
distribution w £ 'D'{Q) if it has a representative {u^)e£{o,i] G £m{^) such that 

u^^w inV'in) ase|0. 

We write U ~ w and call w the associated distribution of U provided U is associated 
with w. 

Regularity theory for linear equations has been based on the subalgebra 
G°°{^) of regular generalized functions in ^(51) introduced in [21]. It is defined 



Wave equation with a discontinuous coefficient depending on time 



5 



by aU elements which have a representative iu'^)ee{o,i] with the property that, for 
all K il, there exists p > such that, for all a £ Ng, 

sup \d^u%x)\ ^ 0{e-P) aseiO. 

Wo observe that all derivatives of u'^ have locally the same order of growth in e > 0, 
unlike elements of fj\/(f2). This subalgebra 0°°{n) has the property Q°°{rt) D 
= C°°{ft), see [21] Theorem 25.2]. Hence, for the purpose of describing 
the regularity of generalized functions, 0°°{il) plays the same role for Q{il) as 
C°°{ft) does in the setting of distributions. The C/°°-singular support (denoted by 
singsuppgoo) of a generalized function is defined as the complement of the largest 
open set on which the generalized function is regular in the above sense. 

We end this section by recalling the notion of slow scale nets. A net (r'^)gg(o,i] 
is called a slow scale net if 

|r^|f = Cl(e-i) aseiO 

for every p > 0. A positive slow scale net is a slow scale net {r'^)ee{o,i] such that 
> for all e £ (0, 1]. We refer to [17] for a detailed discussion of slow scale nets. 

Example 2.3. Let ip he a fixed element of C^(R) such that ip is symmetric, <p' > 
on [—1,0], supp(p C [—1,1] and J^(p{x)dx — 1. Put <Pe(x) — (p(xle)/£. Then 
U £ ^/(M) defined by the class of {fe)e^{o,i\ is associated with the delta function (5, 
and singsuppgcxj U = {0}. On the other hand, if [/ £ Q{M.) is defined as the class of 
(<<5ft,(e))eg(o,i], where (l//i(£))£g(o,i] is a positive slow scale net, then it is associated 
with the delta function again, but singsupp^oo [/ = 0. More generally, for any 
distribution / £ 2?' (17), there exists a generalized function U £ Q°°(yi) which is 
associated with /, see e.g. H]. Thus, any distribution on can be interpreted as 
an element of C/°°(ri) in the sense of association. 



3. Existence and uniqueness of generalized solutions 

We rewrite problem (|l.ip in the form 

d^U - C^dlU = in ^([0, oo) x R), 

f/|t=o = f/o, dtU\t=o^Ui ina(K) ^^'^^ 

in the space of generalized functions, where C is an element of ^([0, oo) x M). In 
the Colombeau setting, existence and uniqueness follows, for example, from results 
in [191 [20], provided the coefficient C is of logarithmic type, i.e., satisfies bounds of 
type 0(log [el) in (|2.ip . When the coefficient depends on time only, this hypothesis 
is not required, as we are going to show in the following existence and uniqueness 
theorem for problem p.ip . 

Theorem 3.1. Assume that C £ f/([0, oo) x M) has a representative (c'^(i))eg(o,i] 
independent of x and satisfying the following two conditions: 
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(i) there exist two constants cq, ci > such that, for any e £ (0, 1] and t > 0, 

ci > c'(f) > Co > 0; 

(ii) for any e G (0, 1], 

\{c^Y{t)\dt < oo. 



Then for any initial data Uq, Ui G ^/(R), problem p.ip has a unique solution 
U e e([0,oo) X R). 



Proof Put V ^dtU - Cd^U and W = dtU + Cd^U. Then problem can be 
rewritten as the Cauchy problem for a first-order hyperbolic system 

{dt + Cd,j:)V = MV - MW in ^([0, oo) x M), 

{dt - Cd^)W = MW - MV in g([0, oo) x M), 

V\t^o = = t/i - {C\t=o)Ui in ^^"^ ^ ■ > 



W\t=o = Wo = Ui + {C\t=oM 



m 



where M = C"/(2C) g 0([O,oo) x M). If problem has a unique solution 

{V, W), then so does problem p.ip . To prove the existence of a solution (y, W), 
let (t;^, w^) be the unique C°°-solution to the Cauchy problem 



[dt + c'{t)d^)v'' = /i=(t)u- - ii''{t)w^, t>Q, .T e 

(a* - c'^(t)9:r)u;'' = n%t)w^ - ^i'' {t)v'' , t>0,xe 

v^lt^o = = -c-(0)(ug)', xGR, 

-;«^|t=o = +c^(0)(ug)', xgR, 



(3.3) 



where (wo)ee(o,i]' ('"i)ee(o,i]' (c^)ee(o,i] and (Ai^)ee(o,i] are representatives of Uq, 
U\, C and M, respectively, such that c*^ is as in the statement and /i'^ = (c^)'/ (2c^). 
For the existence of such {v'^ , w'^), see [20]. If we show that ('y^)ee(o,i] and (w'^)gg(o,i] 
belong to £m{[0,oo) x R), their equivalence classes in ^([0,oo) x R) will form a 
solution of problem p.2p . To show that the zeroth derivatives of v'^ and w'^ satisfy 
estimate (|2.ip . consider the characteristic curves j^{t, x, r) and 7l(i, x, r) passing 
through {t, x) at time t = t which satisfy 

dr'y+{t,X,T) = C^t), '-fl{t,X,t) = X, 

drl^_{t,x,T) = -c^{t), 'yl{t,x,t) ^ X. 
Along these characteristic curves, v' and are respectively calculated as 

w^(i,x) =«g(7;(t,x,0))+ / ^l'{s)iv' ~wn{s,-fl{t,x,s))ds, (3.4) 

Jo 

w^t,x) =w'o{-fLit,x,0))+ I n'{s){w' -v'){s,Y-{t,x,s))ds. (3.5) 
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For each T > 0, we define Kt as the trapezoidal region with corners (0,— 
(T, + ciT), (T, i - ciT), (0, 0- Using and dS^]), we see that 




We add these two incquahties and apply Gronwall's inequality to get 

IK'IU-(Kt) + W\\l^(Kt) < (IkolU-(ifo) + lkolU-(ifo))exp (^2 \n^{s)\dsj . 

On the right-hand side, the terms involving Vq and Wq are of order 0{e^^) for some 
p > 0. The exponential term is uniformly bounded in e by the condition (ii) on . 
Hence, the zeroth derivatives of and w'^ satisfy estimate (|2.ip on Kt- We can 
obtain analogous estimates for all derivatives of and by differentiating the 
equations and using the same argument. Thus, (w^)eg(o,i] and iw'^)ee{o,i] belong 
to £:m([0,cx)) X R). 

For the proof of uniqueness, we only need to obtain the zero-order estimates 
(by Lemma 1.2.3 in [H]), which follow along the same line as above. □ 

We remark that condition (ii) in Theorem 13.11 can be weakened to the re- 
quirement that exp (Jp°° |m'^(^)I dt) = 0{e~P) as e J, for some p > 0. 

4. Propagation of singularities 

In this section we study the phenomenon of propagation of singularities in the 
generalized solution to problem (jl.ip with c{t) — cq -{- (ci — co)H{t — 1), uq = 
and Ml = S, where Co,Ci > 0, Cq 7^ Ci, H is the Heaviside function and S is the 
delta function. Let us begin by regularizing the coefficient and initial data. Let ip 
be as in Example 1 2 . 3 1 and put 

c%t) = (c * ip,){t) = Co + (ci - Co) / H{t-1~ es)ip{s) ds. 

Then c^ c in X''(M) as e J, 0, and the family (c^)£g(o,i] belongs to £m{[0, 00) xR). 
We define C G ^([0, 00) x R) as the class of (c^)eg(o,i]. Then singsuppgcx, C = {t = 
1}. We also define Uq = and Ui G G{M.) as the class of {(pe)e£{o,i]7 where ip again 
is a mollifier as in Example l2.3l (Actually, the moUifier need not be the same as the 
one chosen for the regularization of the coefficient c.) Thus we interpret problem 
(|l.ip with c(i) = Co + {ci — co)H {t — 1) , uq = and ui = 5 as problem p.ip with C, 
Uo and Ui defined above. The existence and uniqueness of a generalized solution 
is then ensured by Theorem 13.11 As may be seen from the following theorem, the 
splitting of the singularities occurs at points of discontinuity of the coefficient (see 
Figure [T]). 
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Figure 1. The t/°°-singular support of the solution U 



Theorem 4.1. Let C, Uq and Ui he as above and let U G Q{[0,oo) x R) be the 
solution to problem ()3.1[) . Then it holds that 

singsuppg^o [/ = .t) I a: = ± /J c(s)(is, t > o| 

U [{t, x) I X = ± (2 c{s) ds - /J c(,s) ds) , t > 1} . ^^'^^ 

Proof. For the sake of presentation, we first assume that ci > cq. To show that 
assertion (|4.ip holds, we first calculate the C/°°-singular support of the solution 
(y, W) to problem (|3.2|) with Vq = Ui and Wq = 0. As may be seen from the proof 
of Theorem 13.11 the solution (V, ly) has a representative (w^, w'^)£g(o,i] i which 
satisfies the integral equations p.4p with Vq = (p^ and (j3.5p with Wq = 0: 

= .^e(7^(t,x-,0)) + / Ai=(s)(u'' - u;=)(s,7^(t,a;,s))ds, 
Jo 

u;'^(t,a-)=: / fi'^{s){w'^ ~ v'^){s,'yt{t,x,s))ds, 
Jo 

where 

7^(i,x,t) = X, 
7I (i, cc, <) = X. 



jlit,X,s)=jlit,X,0)+ C%T)dT, 

Jo 

7i(t,x,s) =7i(t,x,0)- / c^(r)rfr, 



Wave equation with a discontinuous coefficient depending on time 



9 



The solution is obtained by iteration, see [21]. From this, we can check 

that v^{t, x) > and w^{t, < for t > and x G M. Using this fact, we get 

v'it,x)>ipehl{t,x,0))>0, (4.2) 

w'^{t,x)<- / fi'^{s)v'^{s,jt{t,x,s))ds 
Jo 

< - [ f/is)if,iY+{s,r-{t,x,s),0))ds 
Jo 

= -J^ tJ-''{s)(p, (^r_{t,x,0) -2 c^{T)dT^ ds <0. (4.3) 
Step 1. We first prove that 

singsuppgoo F = |(t, a;) I X = J c{s)ds,t>0^. (4.4) 



It is easy to check that V g ^([0, oo) xK) vanishes off {(i, x) \ x = c(s) ds, t > 0} 
and thus is ^""-regular there. We will show first that {{t,x) \ x = /J c(s) ds, t > 0} 
is contained in singsupp^oo V. 

Fix t > arbitrarily, and consider 



At 



v'{t,Y+iO,e,t)) - v'{t,Y+i^,0,t)) 



Y+{o,s,t)-T+io,o,t) 



By inequality (|4?2| . we have {t,jl{0,0,t)) > (p{0)/e > 0. Noting that v^{t,x) ^ 
for a; > ■y^{0,e,t), we get w^(t,7^(0,e,t)) = 0. Furthermore, 7^(0, e,t) — 
7^(0, 0,i) ~ e. Hence, we see that At^e > '■Pi^)!^'^- The mean value theorem shows 
that there exists x\ € (7^(0, 0, t), 7^(0, e, t)) such that 

\d.v^{t,x\)\ > 

Repeat this process to find that for any n>2 there exists xf^ € {x'^_i,j'^{0,e,t)) 
such that 

|a>-«,<)l>|i. 

Since — > c{s) as e J, 0, assertion (|4.4p holds. 
Step 2. We next prove that 



singsuppgcx, = |(t,.T) I 2; = 2y c{s) ds — J c{s)ds,t>l^. 



(4.5) 



The proof is similar. We can easily check that W equals outside {{t,x) \ x = 
2 c(s) ds ~ c(s) ds, t > 1} and thus is C7°°-regular there. We will show below 
that {{t, x) \ X ~ 2 c(s) ds — J* c(s) ds, t > 1} C singsupp^oo W. 
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Fix t > 1 arbitrarily. Let e < t — 1 , and consider 

w'it,YMJo c'jr) dT,t)) - w'{t,Y_{l -£,-£ + /o ''c^t) dT,t)) 



7l (1, /; (r) dT, t) - 7i (1 - £, -£ + /o ' c-(t) dr, 
We see that, for s > 0, 



7l(^l,^ c'{T)dT,S^ ^2 c'{T)dT- C%T)dT. (4.6) 

Using this, inequahty (|4.3p and supp^^ C [1 — £, 1 + e], we get 

(t,jt (l, rc'{T)dT,t)] < - C^\'{S)^, (2 [\%T)dT) ds. (4.7) 



Choose a £ (0,2ci) so that ip{a) > 0. Then, for 1 - e < 1 - (ae)/(2ci) < s < 
1 + (ae)/(2ci) < 1 + e, we have |2 c^(t) dr] < ae. Using this and noting that 
ip' > on [—1, 0] and the symmetry of (ys, we find that (y5e(2 c'(r) dr) > ip^{ae) 
for 1 — (ae)/(2ci) < s < 1 + (ae)/(2ci). Therefore, by inequahty (j4.7p . we obtain 

w;M i,7i 1, / c^(T)dT,t < - / fi%s)ip,{ae) ds 

\ \ Jo J J Jl~(ae)/{2ci) 

^\ c^l + (a£)/(2ci)) ^{a) 
2 °Sc-(l-(a£)/(2ci)) ■ e ' 

where 

^ c^l + (a£)/(2ci)) CO + (ci - cq) J!^''^' ^{y) dy 



c^(l - {as) /{2c,)) + _ co)/_-:/^'^^^ 9^(2/) dy 
is independent of e and greater than 1. Thus, we obtain 

We also see that, for s > 0, 

7l^l-£,-e + ^ c%T)dT,s^ = -e + 2 d'{T)dT ~ d'{T)dT, (4.8) 
on which = 0. Furthermore, by ()4.6|) and ()4.8p . we have 

7! fl, / c''{T)dT,t^ - (^l - e,-e + J c%T)dT,t 

cf{T)dT < (1 + 2ci)e. 

Hence, -64. ^ > (log /3-(y9(a))/(2(l + 2ci)£^). By the mean value theorem, there exists 
xf £ (71(1 —£,—£ + Jq " c^{t) dr, t), 7! (1, /q^ (^^{t) dr, t)) such that 
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Repeating this process gives that for any n > 2, there exists a;^ G (7l(l —£,—£ + 
Jq ^ c^(t) dr, t),x^^_i) such that 

Since — >■ 2 c(s) ds — c(s) as e J, 0, assertion (|4.5p holds. 

Step 3. We finally prove that assertion (|4.ip holds. Similarly to Steps 1 and 2, we 



can show that 

'' ^ c(s) ds ) , t > 1 }> , (4.9) 



sing suppgoc y = I (i, x) I X = — ^2 J c(s) ds — J i 
singsuppgoo = |(i, x) I a- = — y c(s) ds, t > ol , 



(4.10) 



if Vb = and Wq ~ Ui. By the definitions of V and W, the C/°°-singular support 
of U coincides with the union of (|0)) . ((43)) . (|49|l and (|4?T0)) . Thus, assertion (|4?T|) 
follows. 

The case cq > ci can be treated by the same arguments, using a change of 
sign in J7i and invoking the linearity of the equation. The proof of Theorem 14. II is 
now complete. □ 

Wc remark that assertion (|4.ip in Theorem 14.11 still holds when Ui € tJ(M) 
is given by the class of ((/3")£g(o^i] with n £ N, i.e., Ui is any power of the delta 
function. 



5. Slow scale coefficients and regularity along the refracted ray 

We here discuss how the regularity of the coefficient C affects that of the solution 
U to problem (j3.ip . As mentioned in ExamDlc l2.3[ one may regularize the pieccwise 
constant propagation speed c{t) ~ cq + (ci — CQ)H{t — 1) in such a way that the 
corresponding element C € S(K) belongs to tJ°°(R). In fact, it suffices to use a 
moUifier (ph{e)i where (l//i(e))£g(o,i] is a positive slow scale net. It has been shown 
in [18| that a bounded element C of Q{M.) belongs to Q°°{M.) if and only if all 
derivatives satisfy slow scale bounds. Consider problem p.ip with Uq = and Ui 
given by the class of (<^e)ee(o,i] as in Theorem 14. II Then as may be seen from the 
following theorem, the refracted rays 

{t,x) \ X = 2 / c(s) ds - c(s) ds, t > 1 
Jo Jo 
= {{t,x)\x = co~ci{t-l), t>l} 

r+ = ((t,a;) I X = -2 / c{s)ds+ [ c{s)ds,t>l 



= {it,x)\x^-CQ + ci{t~l), t>l} 
do not belong to singsuppgoc U, if C belongs to tJ°°(R) (sec Figure [5]). 
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Figure 2. The fj°°-singular support of the solution U with slow 
scale coefBcient 



Theorem 5.1. Let C, Uo and Ui he as above and let U E Q{[0,oo) x R) be the 
solution to problem If C £ ^°°(K), then 

singsuppgoo U ^ l^{t,x) \ X = ± J c{s)ds, t > j- . 

Proof. As can be seen from the proof of Theorem l4.1[ it suffices to consider problem 
([3:2)) when Vo = Ui and Wo = 0: 

{dt + Cd^)V ^ MV - MW in a([0,oo) X R), 

{dt ~ Cd.^)W ^ MW - MV in a([0,oo) X R), (5.1) 

V\t=o = Ui, W\t=o = in g{R). 

The same argument as in the proof of Theorem 14.11 shows that 

c{s)ds, t > 



sing suppgoo F = I (t, a;) \ J 



holds. Thus, it remains to show that W is C/°°-regular in a neighborhood of r_. 

A simple calculation shows that the commutator oi dt~C{t)dx and dt+C {t)dx 
is given by 

[dt - C{t)dx,dt + C{t)dx] = 2C'{t)dx. 
Using this commutator relation and the second line of system (j5.ip . we have 

{dt - cdx){dt + cdx)w = {dt + cdx){dt - cdx)w + 2C'dxW 

= {dt + Cdx){M{W - V)) + 2G'dxW 

= M{dt + Cdx)W ~ M{dt + Cdx)V + M'{W ~V) + 2C'dxW. 
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From the first line of system (|5.1[) , 

{dt+Cd,)V ^M[V -W). 

Further, 

d^w = ^ ((a* + cd.,)w ~ [dt - C9,)w^) = ^ ((a* + c9,)vf + m{v - w) 

Together with C'/(2C) = M this leads to 

{dt - + cd^)w = 3A/(at + cd,)w + (m^ - M')(y - 1^). 

Let i^T be as in Theorem 14.11 Then the representatives of (V, W^) are of 

order e"-*^ on Kt for some N. Further, M and its derivatives are slow scale by 
assumption; hence Ka*^)^ — (m^)'I < ^li^) for some slow scale net (Ai(£))eg(o,i]. 
From integration along characteristics and Gronwall's inequality as in the proof of 
Theorem 13. 1[ we get 

m + c'd,)w'\\L^^K^^ < Ai(e)£-^Texp 3\^j.%s)\ds^ = C)(Ai(£)e-^) 

since \ fj.'^{s) \ ds is bounded by construction. Recursively we find that 

{dt ~ Cd^Kdt + Cd^fW = {2k + l){dt + Cd^fW + L{V, W, M, C) 

where L depends hnearly on V, W, {dt + Cdx)W, . . . ,{dt + CdxY~^W and may 
contain derivatives of C, M and their products and powers. The linearity of L in 
the first variables and the slow scale property of C, M and their derivatives yield 
by induction that 

\\{dt + c'dxfw'\\L^(K^) = 0(Afe(e)e-^) 

for some slow scale net (A/c(£))ee(o,i]- 

For the remainder of the proof we may use the fact - already shown in 
Theorem 14. II - that V vanishes on r_, hence 

{dt - Cdx)W = MW, {dt - Cdx)^W = M{dt ~ Cdx)W + M'W, . . . 

along r_. Again the higher order directional derivatives of W in the direction 
dt — Cdx can be estimated inductively and bounded by some slow scale net times 
. Similarly, all mixed derivatives {dt ± Cdxf{dt T CdxYW can be estimated 
by reduction to previously computed terms. In conclusion, all derivatives of W 
can be bounded by some slow scale net times £^^, hence are of order 0{e^^^^). 
Thus W is along r„. □ 

We remark that the result can be generalized in various ways. For exam- 
ple, the integral bound on |^^| can be replaced by the requirement that the net 
exp (Jg°° dt) is slow scale. The proof is expected to go through for arbitrary 

initial data Ui whose 5°°-singular support is {0}. 
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6. Associated distributions 

We consider again the wave equation in one space dimension with propagation 
speed c = c{t) depending on time, i.e., 

d^u - c{tfdlu = (6.1) 

with initial conditions 

u\t=o = uo, dtu\t=o = ui. (6.2) 

If c{t) is piecewise constant, say c{t) = cq + (ci — co)H{t — 1), and mq, ui G C', this 
problem has a unique solution 

MeCi([0,oo) : P'(]R)) n (C2((0,l) : P'(M)) ©C2((l,oo) :23'(K))). (6.3) 

Thus the problem is interpreted as a transmission problem across a discontinuity 
at t = 1; the transmission condition is that the solution should be a continuously 
differentiable function of time with values in 2?'(]R). It is simply obtained by solving 
the wave equation for t < 1 and for t > 1 and taking the terminal values aX t ~ 1 
as initial values for t > 1. 

On the other hand, imbedding tio, ui and c into the Colombeau algebra 
by convolution with suitable compactly supported moUifiers (concerning c, we 
use a moUifier as in Example 12. 3p . problem (j6.1l) . (|6.2p has a unique solution 
U S 0([O, oo) X R). We are going to show that the Colombeau solution is associated 
with the piecewise distributional solution. 

Theorem 6.1. Let c{t) he a piecewise constant, strictly positive function and let 
Uq,Ui G 2?'(]R). Then the corresponding generalized solution U G t/([0, oo) x R) is 
associated with the piecewise distributional solution u satisfying li6.S]] . 

Proof. We first assume that wq G C^(]R) and ui G C"'^(R) and that both have com- 
pact support. Fix some T > 1. Let {u'^)e^(a.i\ be a representative of the generalized 
solution that vanishes for x outside some compact set, independently of i, < t < 
T. This and the smoothness imply that belongs to C°°([0,T) : iJ°°(R)). Thus 
we can use energy estimates. Multiplying the wave equation u^^ — c'^(t)^u|^ — 
by Uj and integrating by parts we get 

/ (« + (c")'«t) dx = 0. 
Js. 

Observing that 

~{c^u%f = {c^fulul, + c^{c^)\ulf 

we obtain that 

for t G [0,T) and thus the energy estimate 

/ (Kp + (c-)2|<p)dx< / (|«fp + (c^)2KJ2)dx + 2 fw{c^)'\ ( \ul\^dxdt. 
Jvl Jr Jo Jr 
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By assumption, c^(c^)' is bounded in L^{0,T), and so Gronwall's inequality shows 
that |u|(t, x)p and hence also J^\u^(t, x)\'^ dx remain bounded on [0,T], 
uniformly in e. In particular, the family {u'^)e(£{o,i] is bounded in C([0, T) : iJ^(R)) 
as well as in C"'^([0,r) : i^(R)). By construction, the supports of the functions 
u"^ are contained in a common bounded set. Thus the first property implies that 

t))eg(o_i] is relatively compact in L^(R) for every t. The second property 
shows that (u'^)£g(o,i] is an equicontinuous subset of C([0,T) : L^(]R)). By Ascoh's 
theorem, the net (u'^)£g(o,i] is relatively compact in C([0,T) : L^(E)). 

Taking x-derivatives in equation (|6.ip we can apply the same argument to 
m| and conclude that both a;)p dx and Jjj {u'^^it, x)\'^ dx remain bounded 

on [0,r], uniformly in e. By the differential equation (j6.ip . the same is true of 
\Uff{t,x)\'^ dx. By the same argument as above, (ut)eg(o,i] is relatively compact 
in C{[0,T) : L'^{m.)) and so {u''),^(^o,i] is relatively compact in C^{[0,T) : L'^{R)). 

There exists a subsequence {u'^'')k such that 

lim u^" ^u(^CH[Q,T):L^{M.)) C C^([0, T) ; 2?'(M)). 

fe— s-oo 

But on every compact subinterval of (0, 1) and of (1,T), is identically equal to 
Co or ci, respectively, when e is sufficiently small. This implies that u is a distri- 
butional solution of the wave equation (|6.ip on both strips. Since u G C([0,T) : 
Z)'(R)), so is Uxx- From the equation, we get that 

u G (C2((0, 1) : V'{R)) © C2((l, oo) : V'{m))) . 

In other words, u is the unique piecewise distributional solution to (|6.ip . (|6.2p . 
Consequently, the whole net {u'^)e^{o,i\ converges io u = u. 

If ito G C^(M) and ui £ C^(]R) do not have compact support, we consider an 
arbitrary rectangle [— i?, B] x [0, T] and take a cut-off function x(2;) identically equal 
to 1 on a neighborhood of [—R — rmaxo<t<T c(<), i? -I- Tmaxo<t<T c(t)] x [0, T]. 
By the derivation above, the Colombeau solution with initial data xmq, X^i is 
associated with the piecewise distributional solution with the corresponding initial 
data. But both solutions coincide with the corresponding solutions without cut-off 
on the rectangle \—R, B\ x [0, T], by finite propagation speed. This shows that the 
association result holds without the assumption of compact support. 

In the next to last step, take mo,mi G 2?'(R) and assume that they are dis- 
tributions of finite order. Write uq = Uq + u^^ where Uq has its support bounded 
from the right, and Uq has support bounded from the left, and similarly for ui. 
There is an integer n such that Uq * In £ C^(R) and -uf * /„ G C^(M), where /„ 
is the n-fold convolution of the Heaviside function with itself; a similar assertion 
holds for Uq * In and Uq * /„. We may assume without loss that Mq ^ =0. Let 
u and U be the piecewise distributional solution and the Colombeau solution with 
initial data uq,ui, respectively. Let u and U be the piecewise distributional solu- 
tion and the Colombeau solution with initial data uq * /„ and ui* In, respectively. 
By the previous step, U is associated with u. But 9"C/ = U and d^u = u, since 
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both satisfy the wave equation with initial data uq, ui in the appropriate settings. 
Therefore, U is associated with u as well. 

Finally, if uq,ui £ P'(R) are arbitrary distributions, we use the fact that 
they are locally of finite order and argue by cut-off and finite propagation speed 
as above. □ 

Remark 6.2. The compactness argument in the first step of the proof can be set 
up in various spaces. For example, the uniform boundedness of /jj a;)p da; 

on [0,r] implies the equicontinuity of (w'^)e6(oa] in C{[0,T) : H^{R)) C C([0,T) : 
C(M)). Hence iu^)e(^{o i] is actually relatively compact in C([0,T) : C(M)) = C(R x 

[o,r)). 
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